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This research resolves a long-standing problem on the electromagnetic response of self¬ 
complementary metallic screens with checkerboardlike geometry. Although Babinet’s principle im¬ 
plies that they show a frequency-independent response, this unusual characteristic has not been 
observed yet due to the singularities of the metallic point contacts in the checkerboard geometry. 

We overcome this difficulty by replacing the point contacts with resistive sheets. The proposed 
structure is prepared and characterized by terahertz time-domain spectroscopy. It is experimen¬ 
tally conhrmed that the resistive checkerboard structures exhibit a hat transmission spectrum over 
0.1-1.1 THz. It is also demonstrated that self-complementarity can eliminate even the frequency- 
dependent transmission characteristics of resonant metamaterials. 

BAGS numbers: 78.20.Ci, 78.67.Pt, 42.25.Bs, 64.60.ah 


Duality is one of the key concepts in physics and en¬ 
gineering, as exemplified by electromagnetic duality [1], 
T duality in string theory [2], and the duality of elec¬ 
trical circuits [3]. It relates two seemingly different sys¬ 
tems or quantities, and sometimes helps us to indirectly 
gain physical insight into intractable problems. A sys¬ 
tem is said to be self-dual if it coincides with its own 
dual. Self-duality is symmetry with respect to duality 
transformations. Problems with self-duality often have 
simple analytical solutions due to the constraints im¬ 
posed by their symmetry. These results are universal 
and do not depend on the details of the problems. For 
example, self-dual symmetry has been utilized to deter¬ 
mine critical temperatures of two-dimensional Ising mod¬ 
els (Kramers-Wannier duality) [4]. Moreover, there are 
examples in electromagnetic systems ranging from dc to 
radio frequency: evaluating the effective conductivity of 
two-phase composite media (Keller-Dykhne duality) [5], 
obtaining constant-resistance electrical circuits [6], and 
making broadband antennas [7]. Recently, self-dual sym¬ 
metry has been applied to the design of metamaterials 
with zero backscattering [8]. 

Let us now focus our attention on the duality and self- 
duality of planar structures composed of two elements, 
which can be represented by two-tone patterns. Inter¬ 
change of these two colors is considered to be a duality 
transformation. If patterns are invariant under the inter¬ 
change of colors, as shown in Fig. 1(a), they are self-dual. 
This type of symmetry is referred to as color symmetry [9] 
or self-complementarity, and it is often found in designs 
of traditional garments and in the impressive works of 
the graphic artist M. C. Escher. 

In optics, there is a well-known duality relationship 
called Babinet’s principle [10], which relates the fields 
scattered by a metallic screen with those scattered by its 
complementary screen, which is obtained by interchang¬ 
ing the areas of metal and the holes. With Babinet’s prin¬ 
ciple, it has been shown theoretically and experimentally 


that an antenna with a self-complement ary shape ex¬ 
hibits a frequency-independent input impedance [7, 11]. 
Similarly, it also predicts that self-complementary metal¬ 
lic screens exhibit frequency-independent responses as 
described below. Here, we will consider the scattering 
problems of self-complementary metallic screens such as 
the “ideal” checkerboard geometry shown in Fig. 1(a). 
Babinet’s principle ensures that the power reflectance R 
of the original problem is equal to the power transmit¬ 
tance Tc of the dual problem (Tc = R). On the other 
hand, the power transmittance T of the original problem 
must be equal to Tc, due to the self-duality (T = Tc). 
Then, the energy conservation law (T + T = T + Tc = 1) 
gives the result that the power transmittance of electro¬ 
magnetic waves through the checkerboard screen is equal 
to 1 /2 and is independent of the frequency of the incident 




FIG. 1. (color online), (a) Examples of planar patterns 
with self-complementarity, (b) The ideal metallic checker¬ 
board structure with self-complementarity, and (c) the resis¬ 
tive checkerboard structure. 
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waves [12]. However, this result seems strange for the 
following reasons: (i) the ideal checkerboard structure 
shows a frequency-independent response in spite of the 
fact that its geometry has periodicity or a characteristic 
length; (ii) the frequency-independent spectrum violates 
Foster’s reactance theorem in the long-wavelength limit, 
which states that the reactance of passive and lossless 
systems must strictly increase with frequency [13]. 

During the past few decades, several experimental and 
numerical attempts have been made to observe the un¬ 
usual frequency-independent transmission spectrum of 
the ideal checkerboard structures [14-16], but none of 
them has succeeded in observing the flat spectrum pre¬ 
dicted by this theory. This inconsistency is attributed to 
the point contacts at the corners of the metal patches in 
the ideal checkerboard geometry illustrated in Fig. 1(b). 
In principle, it is impossible to realize such ideal point 
contacts. Thus, the corners of actual structures have to 
be either connected or disconnected. It is known that 
such metallic objects that are nearly touching exhibit 
singular electromagnetic responses [17]. In addition, as is 
the case of the dc electrical conduction of the system [18], 
the connectivity of the metal corners has a significant in¬ 
fluence over the scattering characteristics at higher fre¬ 
quencies [16, 19, 20]. This critical behavior has been ex¬ 
plained by percolation theory with the identification of 
the ideal checkerboard as the structure representing the 
percolation threshold between connected structures and 
disconnected ones [12]. In the case of self-complementary 
antennas, the connectivity of the metal is not considered 
as a critical problem, because the point contacts work 
as feed points to connect external circuits. Therefore, 
the existence of a flat transmission spectrum in metal¬ 
lic checkerboard screens is still controversial and worth 
pursuing. 

In this Letter, to resolve the long-standing problem 
concerning self-complementarity, we propose using yet 
another intermediate state between connection and dis¬ 
connection instead of the singular point contacts. To be 
more specific, we will replace the corners of the metal¬ 
lic checkerboard structure with resistive sheets, as shown 
in Fig. 1(c). By controlling the resistance of the sheets, 
we can realize the intermediates between the connected 
states and the disconnected ones [21]. We then intro¬ 
duce an extension of Babinet’s principle to deal with the 
resistive elements and show a self-dual condition on the 
resistance. Finally, we experimentally demonstrate in the 
terahertz regime that the checkerboard structure loaded 
with resistive elements exhibits the predicted frequency- 
independent response under this specific condition on the 
resistance. 

The familiar form of Babinet’s principle relates the 
scattering due to a thin metallic structure to that due 
to its complement [10]. Note that the structures are as¬ 
sumed to be perfect electric conductors. To discuss the 
effect of the resistive elements in the resistive checker¬ 


board structure, we need to extend Babinet’s principle 
to finite sheet impedance. Let us consider the problem of 
plane-wave scattering by a structure placed in a vacuum 
for which the spatial distribution of the sheet impedance 
Z{x,y) is as shown in Fig. 2(a), where (x, y) gives the co¬ 
ordinates on the structure. Next, we construct the dual 
problem, shown in Fig. 2(b), where the complementary 
sheet-impedance distribution Zc{x,y) is defined as 


Zc{x,y) 


{Zo/2f ^ 

Z{x,y) ' 


( 1 ) 


Here, Zq ^ 377 Q is the impedance of a vacuum. In 
simple terms, the sheet impedance is inverted at each 
point on the structure. We note that the polarizations of 
the incident waves of the two problems differ by 90 deg. 
Babinet’s principle for finite sheet impedance relates the 
transmittance of these two situations [22, 23]: 


i{uj) + tc(cj) = 1, 


( 2 ) 


where i{uj) and ic{uj) are the complex amplitude trans¬ 
mittance of the zeroth-order diffraction mode in the orig¬ 
inal problem and its dual, respectively, and uj is the an¬ 
gular frequency of the incident wave. Here, the zeroth- 
order diffraction mode refers to the mode that has the 
same wave vector and polarization as the incident one. 
It is easy to confirm that the extended version of Babi¬ 
net’s principle includes the conventional one by recalling 
that Z = 0 and Z = oo correspond to perfect electric 
conductors and holes, respectively. We emphasize that 
Babinet’s principle relates the amplitude coefficients, not 
the power transmittance. In cases with neither polariza¬ 
tion conversion nor diffraction, the power transmittance 
is related (T + Tc = 1) by a combination of Babinet’s 
principle (Tc = R) and energy conservation {T R= 1). 

We apply this principle to the resistive checkerboard 
structure. If the sheet impedance of the resistive sheet 



FIG. 2. (color online). Two problems of plane-wave scat¬ 
tering related through Babinet’s principle extended to finite 
sheet impedance, (a) The original problem, (b) The dual 
problem with the complementary sheet-impedance distribu¬ 
tion and the incident polarization orthogonal to the original 
one. The symbols E, H, and k represent the electric field, 
the magnetic field, and the wave vector, respectively. 
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is equal to Zo/2, the structure is self-complementary, be¬ 
cause resistive sheets are invariant under the transfor¬ 
mation in Eq. (1) and the complementary structure can 
overlap the original one. Hence, we see that the dual 
problem of the resistive checkerboard structure is iden¬ 
tical to the original one when linearly polarized plane 
waves are normally incident. Consequently, we obtain 
i{uj) = tc(co’), and combining this with Eq. (2) leads to 

f(w) = fc(w) = (3) 

This indicates that the resistive checkerboard structure 
shows a frequency-independent transmission spectrum if 
the sheet impedance of the resistive sheets is equal to 
Zo/2. The same applies to the reflection spectrum, and 
we obtain an amplitude reflectance r{uj) = —1/2. There¬ 
fore, half the incident power must be diffracted or ab¬ 
sorbed. In particular, in the long-wavelength limit where 
there is no diffraction, half the power is absorbed by the 
resistive sheets. The theoretical general sufficient condi¬ 
tions for the frequency-independent response have been 
presented in our previous paper (see Ref. [23]). It should 
be noted that Zo/2 is replaced with Zo/(2n) when the 
structure is surrounded by an isotropic dielectric medium 
with refractive index n. 

We will now explain our experimental demonstration of 
the above situation. The resistive checkerboard structure 
was fabricated on a c-cut sapphire substrate (20 mm x 
20 mm X 900 /im, Kyocera) by using the standard pho- 



connected 


disconnected 


FIG. 3. (color online). (a) The unit cell of the resis¬ 
tive checkerboard structure. The dimensions are as follows: 
a = 150/xm, d = 50/xm, w = 30/xm, and I = 30/xm. The 
gray dotted lines indicate the area of overlap of the two lay¬ 
ers. (b)-(d) Photomicrographs of the fabricated checkerboard 
structures. 


tolithography and lift-off technique. Eor the case of nor¬ 
mal incidence, the refractive index of the c-cut sapphire 
plate is that for ordinary waves, i.e., nsa ^ 3.1 in the 
terahertz regime [24]. The unit cell of the fabricated re¬ 
sistive checkerboard structure is shown in Eig. 3(a), and 
its photomicrograph is shown in Eig. 3(b). The struc¬ 
ture consists of two layers. A resistive titanium (Ti) 
layer 19 nm thick and a conductive aluminum (Al) layer 
400 nm thick were deposited by electron-beam evapo¬ 
ration at room temperature. The Ti rectangles have 
a margin of d — I = 20 /xm that overlaps with the Al 
patches to assure electrical contact. With the terahertz 
time-domain spectroscopy (THz-TDS) technique (see the 
Supplemental Material for the experimental details [25]) 
and Tinkham’s equation [27], the sheet impedance of the 
19-nm-thick Ti film was estimated to be 0.98 xZq/ (2nsa) • 
The thickness of the Al layer was determined by taking 
into account the skin depth of terahertz radiation in Al, 
^100 nm at 1 THz [28]. Eor comparison, the connected 
and disconnected checkerboard structures were also fab¬ 
ricated as shown in Eigs. 3(c) and 3(d). Their dimensions 
are the same as the resistive one, although they comprise 
only an Al layer and do not have self-complementarity. 

To investigate their transmission properties, the fab¬ 
ricated structures were also characterized with THz- 
TDS [25]. The probe terahertz beam was linearly polar¬ 
ized in the y direction and focused on the samples under 
normal incidence with respect to the sample plane. In or¬ 
der to apply Babinet’s principle, it is necessary to satisfy 
the mirror-symmetry condition with respect to the sam¬ 
ple plane. During the measurements, the surfaces of the 
checkerboard structures were covered with another plain 
sapphire plate so that the structures were symmetrically 
sandwiched between sapphire plates. We used a pair of 
plain sapphire plates as a reference. The amplitude trans¬ 
mission coefficients i{uj) of the samples were calculated by 

t(cj) E/sample (^)/-^sapphire (^) 5 whcrC 

represents the Eourier transform of a recorded electric 
field of a terahertz pulse transmitted through a sam¬ 
ple (sapphire). In the calculation, the echo pulses caused 
by reflections at the boundaries of the substrate were re¬ 
moved from the temporal waveforms by multiplying by 
a time window. We note that the effect of the diffracted 
waves is negligible if the detector is sufficiently distant 
from the samples, and the terahertz beam can be ap¬ 
proximated by a plane wave in the neighborhood of the 
focal point. Therefore, we can regard the experimental 
results as the theoretical transmission coefficients of the 
zeroth-order diffraction mode. 

The measured amplitude transmission spectra of the 
checkerboard structures are shown in Eig. 4(a) (see 
Eig. S2 in the Supplemental Material for comparison 
with numerical simulations [25]). It is clearly confirmed 
from the result that the resistive case shows a nearly 
frequency-independent spectrum in both the real and 
imaginary parts. Note that the flatness continues beyond 




















4 


(a) 



(b) 


connected 



FIG. 4. (color online). The results of the transmission measurements, (a) The transmission spectra, (b) The loci of the 
amplitude transmission coefficients in the frequency range of 0.1-0.55 THz. 


the diffraction frequency (0.64 THz) or the homogeniza¬ 
tion limit. On the other hand, the spectra of the con¬ 
nected and disconnected checkerboard structures highly 
depend on the frequencies of the incident waves. Fig¬ 
ure 4(b) shows the loci of the amplitude transmission co¬ 
efficients in the frequency range of 0.1-0.55 THz, where 
diffraction is negligible. We can easily confirm that the 
locus of the resistive case stays in close proximity to 
(Ret,Imt) = (0.5,0). On the other hand, as the fre¬ 
quency increases, the loci of the connected and discon¬ 
nected cases move in a clockwise direction along a circle 
centered at (0.5, 0) with radius 0.5. The constraint of the 
motion to the path defined by this circle is due to energy 
conservation in the nondiffraction regime [29]. 

Equation (3) holds for any self-complementarity struc¬ 
tures with n-fold rotational symmetry (n > 3) [23]. 
The consequent frequency-independent response is coun¬ 
terintuitive, especially when these are resonant struc¬ 
tures, which are often building blocks of metamateri¬ 
als. Here, we consider resistive self-complementary struc¬ 
tures loaded with electric-inductive-capacitive (EEC) res¬ 
onators [30], and their complementary structures [31]. 
Their shape and dimensions are shown in Eig. 5(a). They 
are designed to exhibit an inductive-capacitive resonance 
at 0.41 THz when surrounded by sapphire. A photomi¬ 
crograph of the fabricated planar metamaterial is shown 
in Eig. 5(b). The details of the fabrication are the same 
as above. The EEC resonators and their complements 
are placed at the centers of the metallic patches or holes 
in the resistive checkerboard structure to maintain self¬ 
complementarity. The connected and disconnected struc¬ 
tures were also prepared for comparison, as shown in 
Eig. 5(c). We measured the electric fields transmitted 
through the samples by THz-TDS. Eigure 5(d) shows 
the measured waveforms of the electric fields after pass¬ 
ing through the reference sapphire plates and the meta¬ 
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FIG. 5. (color online), (a) The dimensions of the EEG 
resonators and the complementary ones, which are built in 
resistive checkerboard structures: p = 65 //m, s = 8 pm, 
and g — S pm. (b) A photomicrograph of the resistive 
self-complementary metamaterial, (c) Photomicrographs of 
the connected and disconnected checkerboard metamaterials, 
(d) The recorded temporal waveforms of the terahertz elec¬ 
tric fields that were transmitted through a sapphire reference 
and the disconnected, connected, and resistive metamateri¬ 
als. The curves are vertically offset by 2 units for the sake of 
clarity. 


materials (see Fig. S7 in the Supplemental Material for 
their transmission spectra [25]). In both the connected 
and the disconnected cases, the waveforms are highly dis¬ 
torted as compared to those of the reference; this is due 
to the frequency-dependent transmission characteristics 
of the metamaterials. In addition, we can see persis¬ 
tent oscillations after the main pulse, and these are the 
evidence that the incident energy is temporarily stored 
in the resonant structures and then released afterward. 
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On the other hand, the distortion in the resistive case 
is obviously small, and persistent oscillations are not ob¬ 
served. This clearly indicates that self-complementarity 
suppresses the frequency dependence of the transmission 
induced by the resonant structures. We note that resis¬ 
tive loading without self-complementarity does not cause 
the frequency independence (see Fig. S5 in the Supple¬ 
mental Material [25]). 

In conclusion, we performed an experiment that 
demonstrated the controversial frequency-independent 
transmission spectra of self-complementary screens by 
replacing singular point contacts with resistive sheets. 
Consequently, this work revealed that the controversy 
arose from an implicit assumption about the power trans¬ 
mittance, T -f- Tc = 1, which does not hold in self-dual 
cases. We also showed that self-dual symmetry can sup¬ 
press even the strongly frequency-dependent response of 
the resonant structure. The phenomenon reported here 
is universal, because it is based on the self-duality of the 
system. Thus, it is expected to play an important role in 
a broad range of the electromagnetic spectrum and to be 
helpful for practical applications such as designing broad¬ 
band metamaterials. In addition, the flat transmission 
and reflection characteristics of the self-complementary 
screens are applicable to broadband coherent perfect ab¬ 
sorption in artificial structures [23, 32]. We note that 
there is no theoretical limitation on the size of the re¬ 
sistive sheets and that the energy of the incoming radi¬ 
ation can be concentrated on highly subwavelength re¬ 
gions; thus this has the potential to enhance nonlinear 
phenomena and light-harvesting applications. Finally, 
we point out that similar nearly flat transmission spec¬ 
tra in finite frequency ranges have been observed in thin 
metal films close to the metal-insulator transition, where 
metal islands are randomly connected by lossy narrow 
necks [33, 34]. Although these involve imperfections such 
as randomness, broken mirror symmetry, and material 
dispersion, they can be regarded as statistically and ap¬ 
proximately self-dual screens. Thus, the perspective of 
this Letter provides another physical insight into their 
behavior, beyond the conventional view from effective 
medium theory [35]. 
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In this Supplemental Material, we present the results of several numerical studies which support the conclusions of 
the main article. Furthermore, we provide supplemental information on the experiments described in the main article. 
The figures of this Supplemental Material are referred to as “Fig. Sx”. 


I. NUMERICAL SIMULATIONS 
A. Comparison between experiment and simulation 

To estimate the effects of experimental imperfections on the transmission spectra and compare with the experimental 
results shown in Fig. 4, we performed numerical simulations by a commercial finite-element method solver (COMSOL 
multiphysics). The model used in this simulation is shown in Fig. SI. The unit cell is filled with a lossless 
dielectric material whose anisotropic dielectric permittivity tensor e = diag[(nsa )^5 (^Sa)^, (^Sa)^]’ where nsa = 3.1 
and = 3.4 [1]. Strictly speaking, Babinet’s principle does not hold for anisotropic media such as sapphire. 
However, for the case of normal incidence, the c-cut sapphire substrate can be considered approximately as an 
isotropic medium. Metal is assumed to be infinitely thin perfect electric conductor (PEC) for simplicity, and resistive 
sheets are expressed as sheet impedance Zg. It is inevitable in the “sandwich” measurement that the surface of the 
checkerboard is incompletely covered by the dielectric plate. There should be an air gap between the checkerboard 
and the cover plate. To incorporate this experimental imperfection into the simulation model, a thin air layer with 
thickness t = 2 jam is inserted between the cover plate and the checkerboard screen. 

The calculated transmission spectra are shown in Fig. S2(b) for Zg = {10“^, 1,10^} x Zo/{2nsa), which correspond 
to the connected, resistive, and disconnected cases, respectively. The simulations agree with the experimental results 
qualitatively well in spite of the simplified model. 


Top view 

y 



Side view 




FIG. SI. A schematic picture of the unit cell used in the numerical simulations (a = 150/xm, h — 450/xm, d = 30/xm, and 
t — 2 jim ). 
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(a) Experiment 



(b) Simulation 



FIG. S2. (a) The measured transmission spectra of the checkerboard screens, (b) The calculated transmission spectra. 
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B. Effect of sheet impedance on transmission spectra 

Here, we see how the transmission spectrum changes when is varied around the self-dual condition — 
Zo/(2nsa)- Figure S3 shows the model used in the simulations. For simplicity, the dielectric environment is as¬ 
sumed to be isotropic, and the air gap is not considered. The simulations were performed by Ansys hfss. The 
calculated transmission spectra for Z^ = {4“^,2“^, 1,2,4} x Zo/(2nsa) are shown in Fig. S4. From the results, it is 
confirmed that the frequency dependence of the transmission is suppressed as Z^ approaches the self-dual condition. 
Although there is a slight deviation from I = 0.5 even if Z^ = Zo/(2nsa), we confirmed that it can be reduced by 
increasing the mesh density of the model. 



FIG. S3. The model used in the numerical simulations of the checkerboard with the resonant structures (a = 150/xm, h = 
300 /xm, d = 30 /xm, p = 65 /xm, s = 8 /xm, and p = 8 /xm). 



FIG. S4. The calculated transmission spectra of the checkerboard structures with the ELG structures for several values of the 
sheet impedance. 
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C. Resistive structures without self-complementarity 


In order to demonstrate the crucial role of self-complementarity, we calculated the transmission spectra of the 
resistive checkerboard structures without the ELC or complementary EEC (CELC) resonators shown in Eigs. S5(a) 
and S5(b). They are not self-complementary and contain resistive sheets. The model used in this simulation is the 
same as Eig. S3 except the geometry, and is set to Zo/(2nsa)- Eigure S5(c) shows the calculated transmission 
spectra. We can confirm that the structures without self-complementarity do not show the frequency-independent 
response although they contain resistive sheets. This indicates that self-complementarity of the whole structure plays 
an essential role in realizing the frequency-independent transmission spectrum. 


(a) 



without ELC resonators 


(b) 




FIG. S5. (a) The resistive checkerboard structure without the ELC resonators, (b) The resistive checkerboard structure without 
the CELC resonators, (c) The calculated transmission spectra of the resistive checkerboard structures with and without self¬ 
complementarity. 
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II. SUPPLEMENTAL INFORMATION ON EXPERIMENTS 


A. Experimental setup of THz-TDS 


We used the terahertz time-domain spectroscopy (THz-TDS) technique [1] for characterization of the fabricated 
samples. The experimental setup is as follows. Figure S6 shows a schematic picture of the THz-TDS system. 
Photoconductive THz antennas (PC As) on low-temperature-grown GaAs equipped with silicon hyper hemispherical 
lenses (Batop) were used as the emitter and detector. These antennas were excited by a femtosecond fiber laser 
(F-lOO, Imra) with a wavelength of 810 nm, through 1.5m-long optical fibers. The group delay dispersion (GDD) 
introduced by the optical fibers was compensated by grating pairs [2]. The emitted THz beam was linearly polarized 
and focused by Tefion lenses on the samples under normal incidence with respect to the sample plane. The beam 
diameter at the beam waist is estimated as 3.4 mm at 0.5 THz. The transmitted THz electric fields were coherently 
recorded in the time domain. All the THz components were placed in a box purged with dry air to reduce water 
vapor absorption of the THz beam. 
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FIG. S6. A schematic picture of the THz-TDS system (HWP: half-wave plate, PBS: polarizing beam splitter). 
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B. Transmission spectra of checkerboard screens with resonant structures 

Figure S7 shows the transmission spectra which are obtained by Fourier-transforming the temporal waveforms 
shown in Fig. 5(d). While the connected and disconnected cases exhibit complex spectra due to the interference of 
multiple resonances, the resistive case shows a relatively flat one around i = 0.5 in the observed frequency range. This 
indicates in the frequency domain that self-complementarity suppresses the frequency dependence of the transmission 
induced by the ELC resonators. 



Frequency (THz) 

FIG. S7. The measured transmission spectra of the checkerboard screens loaded with ELC resonators. 
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